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Abstract 

Using 1 .8 fb^ 1 of pp collisions at a center-of-mass energy of 7 TeV recorded by the ATLAS de- 
tector at the Large Hadron Collider, we present measurements of the production cross sections of 
T(1S,2S,3S) mesons. T mesons are reconstructed using the di-muon decay mode. Total produc- 
tion cross sections for p T <70 GeV and in the rapidity interval \y r \ <2.25 are measured to be 
8.01±0.02±0.36±0.31 nb, 2.05±0.01±0.12±0.08 nb, 0.92±0.01 ±0.07±0.04 nb, respectively with 
uncertainties separated into statistical, systematic, and luminosity measurement effects. In addition, 
differential cross section times di-muon branching fractions for T(1S), T(2S), and T(3S) as a function 
of T transverse momentum p T and rapidity are presented. These cross sections are obtained as- 
suming unpolarized production. If the production polarization is fully transverse or longitudinal with 
no azimuthal dependence in the helicity frame the cross section may vary by approximately ±20%. 
If a non-trivial azimuthal dependence is considered, integrated cross sections may be significantly 
enhanced by a factor of two or more. We compare our results to several theoretical models of T 
meson production, finding that none provide an accurate description of our data over the full range of 
T transverse momenta accessible with this dataset. 
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cross sections for pt < 70 GeV and in the rapidity interval \y r \ < 2.25 are measured to be 8.01 ± 
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theoretical models of T meson production, finding that none provide an accurate description of our 
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I. INTRODUCTION 

Since the discovery of the J/tp and T mesons [l], 0| 
the study of heavy quark- antiquark systems has provided 
valuable input for our understanding of Quantum Chro- 
modynamics (QCD). However, despite being one of the 
simplest systems in QCD it has proven difficult histori- 
cally to describe the production properties of these states 
adequately. Results on J /ip and Y hadroproduction and 
polarization [H-[l2| exhibit inconsistencies between mea- 
surements and theoretical predictions fl3| . 

Measurements at the Large Hadron Collider (LHC) 
of differential production spectra of various charmonium 
and bottomonium states together with measurement of 
their spin-alignments, prompt double-quarkonia produc- 
tion, and production of quarkonia in association with 
photons, vector bosons or open heavy-flavor final states 
will allow discrimination between different theoretical ap- 
proaches based on color-singlet corrections [T^ - [l6l. c olor- 
octet terms [TtJ , the fcx-factorization approach [18|| and 
other production models [19], and provide additional in- 
put toward an improved understanding of quarkonium 
hadroproduction. 

Studies of bottomonium production complement con- 
current studies of charmonium systems due to the larger 
mass of the bottom quark compared to the charm quark, 
allowing more dependable theoretical calculations than in 
the charmonium family, which suffer from poor pertur- 
bative convergence [171 ] . Extension of cross-section mea- 
surements to higher meson transverse momenta provides 
valuable input to improvements in the theoretical de- 
scription since in this regime different processes can dom- 
inate and, experimentally, the impact of spin-alignment 
uncertainties are mitigated. 

Taking advantage of the large increase in integrated 
luminosity delivered by the LHC in 2011, this paper up- 
dates a previous measurement [8| and reports the T(1S) 
cross section presented differentially in two intervals of 



absolute Y rapidity [20| , \y T \, and fifty intervals of T 
transverse momentum, p^, and as a px-integrated result 
in forty-five bins of absolute Y rapidity. In addition, new 
measurements of the equivalent T(2S) and T(3S) differ- 
ential spectra and their production ratios relative to the 
T(1S) are reported. 

Y production can proceed directly, where the T me- 
son of interest is produced in the hard interaction, or via 
the production of an excited state which subsequently 
decays. This so-called "feed-down" contribution compli- 
cates the theoretical description of quarkonium produc- 
tion as calculation of color-singlct P-wave and higher or- 
bital angular momentum quarkonium production suffers 
from the presence of infrared divergences 21]. From the 
experimental perspective, separation of direct and feed- 
down contributions is hindered by the small mass split- 
ting between the bottomonium states, which impedes the 
detection of additional decay products that would indi- 
cate indirect production. Contributions from feed-down 
vary between the T(1S), T(2S) and T(3S) states due to 
the changing presence of various kinematically-allowed 
decays and influence the inclusive production rate. Study 
of the Y production ratios as a function of kinematic 
variables as presented here thus provides an indirect but 
precise measure of these feed-down contributions. 

This new analysis extends the px range of our previous 
cross-section result Q to 70 GeV, at which new contribu- 
tions to Y production [l5|, [l6| such as from associated 
Y + bb production may play a more important role and 
the impact of the dependence of the production cross 
section on the spin-alignment of the Y is relatively small. 
Both the fiducial cross section, measured in the kinematic 
region with muon transverse momentum p T > 4 GeV and 
muon pseudorapidity |?7 M | < 2.3, and the corrected cross 
section, which is defined in this paper as the cross sec- 
tion in the pT~r] phase space of the Y corrected for the 
acceptance of the decay products to the full Y produc- 
tion phase space, are reported. In the former case, the 



results have no dependence on assumptions about T spin 
alignment, while the latter measurements are more eas- 
ily compared to model predictions and to the results of 
other experiments. 



II. THE ATLAS DETECTOR 

The ATLAS detector [22] is composed of an inner 
tracking system, calorimeters, and a muon spectrome- 
ter. The inner detector directly surrounds the inter- 
action point and consists of a silicon pixel detector, a 
silicon microstrip detector, and a transition radiation 
tracker, all immersed in a 2T axial magnetic field. It 
covers the pseudorapidity [2(| range \rj\ < 2.5 and is en- 
closed by a calorimeter system containing electromag- 
netic and hadronic sections. Surrounding the calorime- 
ters is the large muon spectrometer built with three 
air-core toroids. This spectrometer is equipped with 
precision detectors (monitored drift tubes and cathode 
strip chambers) that provide precise measurements in the 
bending plane within the pseudorapidity range |t;| < 2.7. 
In addition, resistive plate and thin gap chambers with 
fast response times are used primarily to construct muon 
triggers in the ranges |?7| < 1.05 and 1.05 < |ry| < 2.4 re- 
spectively, but are also used to provide position measure- 
ments in the non-bending plane and to improve pattern 
recognition and track reconstruction. Momentum mea- 
surements in the muon spectrometer are based on track 
segments formed in at least two of the three precision 
chamber planes. 

The ATLAS detector employs a three-level trigger sys- 
tem [23| , which reduces the 20 MHz proton bunch colli- 
sion rate to the several hundred Hz transfer rate to mass 
storage. The Level-1 muon trigger searches for hit coin- 
cidences between different muon trigger detector layers 
inside pre-programmed geometrical windows that bound 
the path of triggered muons of given transverse momen- 
tum and provide a rough estimate of its position within 
the pseudorapidity range \rj\ < 2.4. At Level-1, muon 
candidates are reported in "regions of interest" (Rols). 
Only a single muon can be associated with a given Rol 
of spatial extent A(j> x Arj 0.1 x 0.1. This limitation 
has a small effect on the trigger efficiency for T mesons, 
corrected for in the analysis using a data-driven method 
based on analysis of J/i/j — > fj,^ ' \x~ and T — > fi + jjf~ de- 
cays. The Level-1 trigger is followed in sequence by two 
subsequent higher-level, software-based trigger selection 
stages. Muon candidates reconstructed at these higher 
levels incorporate, with increasing precision, information 
from both the muon spectrometer and the inner detec- 
tor and reach position and momentum resolution close to 
that provided by the offline muon reconstruction. 
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III. DATASET AND EVENT SELECTION 



Data for this study were collected during the 2011 LHC 
proton-proton running period between March and Au- 
gust using a trigger that requires the presence of two 
muon candidates with opposite charges that are subject 
to a fit constraining them to a common vertex while tak- 
ing into account track parameter uncertainties. A very 
loose selection on vertex x 2 which is fully efficient for 
signal candidates, was imposed to ensure proper fit con- 
vergence, as well as the requirement of opposite charge 
and that pt^ > 4GeV and < 2.3. This trigger was 
largely unprescaled and collected data at a rate of ap- 
proximately 20 Hz in this period of data-taking. 

Events are required to contain at least one primary 
vertex candidate that has at least five tracks with p^ > 
0.4 GeV, and at least two muons identified by associating 
candidates found in the muon spectrometer with tracks 
reconstructed in the inner detector [j| [22[ . Multiple scat- 
tering in the calorimeters and toroids of the ATLAS de- 
tector degrades the muon spectrometer resolution for low 
energy particles. As the majority of muons selected for 
this analysis have low momentum, we assign values for 
parameters such as pt and rj to the muons based on track 
fits using inner detector information only. To ensure ac- 
curate inner detector measurements, each muon track 
must contain at least six silicon microstrip detector hits 
and at least one pixel detector hit. Muon candidates 
passing these criteria are required to have p!^ > 4 GeV 
and |ry M | < 2.3 and a successful fit to a common vertex. 



Good spatial matching (AR = ^(A4>) 2 + (Arj) 2 < 
0.01) between the muon candidate in both the offline re- 
construction and the trigger used to select the event is re- 
quired to facilitate data-driven estimates of the di-muon 
trigger efficiency. Furthermore, both muons forming an 
T candidate must be associated to their trigger-level can- 
didates in this manner. In this way, the efficiency of the 
trigger requirements on the di-muon candidate are incor- 
porated into the trigger efficiency correction. All di-muon 
candidates passing these criteria are retained for analysis. 



The distribution of the invariant mass of the 
system is shown in Fig. [1] for the selected di-muon candi- 
dates. As is apparent from the plot, the mass resolution 
(120 MeV) for di-muon candidates detected in the central 
region of the detector ( | | < 1.2) is significantly bet- 
ter than the resolution (214 MeV) for those candidates 
falling into the forward region (1.2 < |y w | < 2.25). 
A total of 3.9 x 10 6 and 2.3 x 10 6 candidates with 
8 < < 11.5 GeV are reconstructed in the central 

and forward regions respectively. 
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FIG. 1: The di-muon invariant mass spectum for events used in this analysis. Separate spectra are shown for those events with 
the di-muon candidate (left) in the central region of the detector (|j/ MM | < 1.2) and (right) in the forward region (1.2 < < 
2.25). Overlaid are individual shapes of the fitted T(nS) signals (shaded regions), background-only fit (dashed curve) and the 
total signal plus background shape (solid curve). The background shape is modelled here by a fourth-order polynomial and the 
signal peaks each modelled by a single Gaussian. Also quoted are the fitted mass resolutions of each of the three signal peaks, 
determined from the fit with a common resolution parameter scaling with invariant mass. 



IV. CROSS SECTION DETERMINATION 

Differential T cross sections are measured according to 
the relation: 

x Br(T -> fi+fi-) = -f— ^ — r (1) 

dp T dy y p p ' JCdtx Ap T x A|y| v ' 

where Br(T — > represents the appropriate branch- 

ing fraction of the T(nS) to di-muons, J Cdt is the in- 
tegrated luminosity, Apr and A|y| are the bin sizes in 
T transverse momentum and rapidity, respectively, and 
Nr is the corrected number of observed T(1S), T(2S), or 
T(3S) mesons in a bin. Corrections are applied to ac- 
count for selection efficiencies, bin migration effects due 
to finite detector resolution, and, in the case of corrected 
cross-section measurements, acceptance. 

Determination of the cross sections proceeds through 
several steps. Firstly, a weight is determined for each 
selected di-muon candidate equal to the inverse of the 
total efficiency for the candidate. Secondly, a fit is per- 
formed to the distribution of weighted events binned in 
di-muon invariant mass, m w , to determine the number 
of T(nS) mesons (with n = 1,2,3) produced in each (pt^ 1 , 
y^p) bin. Thirdly, these values are corrected for (ptT, 
y^) bin migrations. Finally, the differential cross sec- 
tion multiplied by the T branching fraction is 
calculated for each state using the integrated luminosity 
and the pt and y bin widths as in Eq. ([1} . 



The weight, w, for each T candidate includes the 
fraction of produced Y — > p + H~ decays with both 
muons falling into the kinematic region pl^ > 4 GeV and 
\rj^\ < 2.3 (the acceptance, A, is used only in calcu- 
lating corrected cross sections), the probability that a 
candidate falling within the acceptance passes the off- 
line reconstruction requirements (the reconstruction ef- 
ficiency, e roco ), and the probability that a reconstructed 
event passes the trigger selection (the trigger efficiency, 
£trig)- The weights assigned to a given candidate when 
calculating the fiducial (u>fid) and corrected (u>tot) cross 
sections are then given by: 

^fid — (^rcco ' ^trig) 3 
Wtot = {A ■ Ercco ' £trig) ■ 

A. Acceptance 

The kinematic acceptance A(pT,y) is the probability 
that the muons from an T with transverse momentum 
and rapidity y fall into the fiducial volume of the detec- 
tor defined by the > 4 GeV and 1 77^ | < 2.3 selection 
applied to each muon in the di-muon pair. In order to 
calculate the acceptance as a function of the rapidity and 
transverse momentum of each of the T states, taking into 
account possible angular dependences in their decays, we 
use an analytic formula describing the decay of T states 
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in their decay frame [2 
d 2 N 



FLAT ATLAS 



dcosO* 



oc 1 + Xg cos^ 9* + \<f, sin^ 9* cos 20* 
+ AflAsin26'* cos 6*. 



(3) 



The helicity frame is used, where 9* is the polar an- 
gle between the /i + momentum in the T decay frame 
and the direction of the T momentum in the laboratory 
frame. The corresponding azimuthal angle 4>* is the an- 
gle between the quarkonium production plane (defined 
in the quarkonium decay frame by the two momenta of 
the incoming protons in that frame) and the quarkonium 
decay plane in the lab frame. 

The Xi coefficients parameterize the spin-alignment 
state of the T. In some parts of the phase space the ac- 
ceptance may depend on these parameters quite strongly. 
We have identified five extreme cases that lead to the 
largest possible variations of acceptance within the phase 
space of this measurement and define an envelope in 
which the results may vary under any physically-allowed 
spin-alignment assumption: isotropic distribution inde- 
pendent of 9* and (f>* (Xg = X^ = Xg^ = 0, labeled 
FLAT); longitudinal alignment (Xg = — 1, = Xg^ = 0, 
labeled LONG); and three types of transverse alignment 
(X e = +1,A = X 8(t> = 0, labeled T+0; A e = +l,X ct> = 



-1.A* 



0, labeled T++; Xg 



-LAw 



-l,A fl 



0, 



labeled T-\ — ). The central values of our measurements 
are derived under the isotropic production assumption. 
The spread in corrected cross-section results derived un- 
der all five assumptions is used to quantify the full en- 
velope of possible variations of the result due to spin- 
alignment. Some constraints on physically-allowed Xi 
combinations exist [241. l25j . In particular, Xg^ must be 
zero for Xg = ±1 and any A^, and |Ag^| < 0.5 for 
Xg = Xcf, = 0. Acceptance corrections with non-zero 
Xg<p were found to lead to corrections with smaller or 
equal variations from the isotropic scenario than any 
parameter-space scenarios on the Ag-A^ plane. Accep- 
tance weights for each polarization scenario, including 
those with non-zero Xg$ can be found in hepdata [26]. 

Figure [5] presents the two-dimensional acceptance map 
for T(1S), as a function of the T transverse momentum 
and absolute rapidity in the unpolarized (FLAT) accep- 
tance scenario as well as the maximum relative variation 
(ratio of largest to smallest acceptance correction) of the 
acceptance maps across all spin-alignment scenarios. 

Acceptance is significantly decreased at very high ra- 
pidities due to the muon pseudorapidity requirement of 
|t7 m | < 2.3. When T are produced at rest, both muons 
are likely to be reconstructed in the muon spectrome- 
ter on opposite sides of the detector. When the T has 
a sufficient boost, the two muons are more likely to be 
found in the same side of the detector and the probability 
that at least one muon is below the required pt^ thresh- 
old increases, leading to a drop in acceptance. At higher 
Pt the T imparts greater momentum to the two muons, 
the likelihood of a muon being below threshold decreases 
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FIG. 2: (Top) Acceptance as a function of Y pt and rapidity 
for the T(1S) in the unpolarized acceptance scenario; (bot- 
tom) the ratio of largest to smallest acceptance correction 
between extreme spin-alignment scenarios. 



and the acceptance asymptotically approaches 100%. As 
T transverse momentum increases the variation of the 
acceptance corrections between spin- alignment scenarios 
decreases, becoming < 10% at the largest pr ranges stud- 
ied. Similar behaviours are observed in the acceptance 
maps for the T(2S) and T(3S). 

This spin-alignment systematic uncertainty only ap- 
plies to the corrected cross-section measurements. For 
fiducial cross sections, measured in the kinematic region 
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p T > 4GeV and \r]^\ < 2.3, we do not have to correct 
our results to the full phase space of produced T mesons. 
Spin-alignment can thus affect only our estimates of the 
reconstruction and trigger efficiencies. Since we measure 
these in many bins of muon px an d 77, effects from dif- 
fering distributions of events over individual bins due to 
different spin-alignment assumptions are negligible com- 
pared to other sources of systematic uncertainty. 



B. Reconstruction and trigger efficiencies 

The efficiency of the offline reconstruction criteria for 
events with muons within the fiducial region is given by: 

Eroco =£trk(PTl,?7l) X £ trk (PT2 , %) ^ 

x ^(Pti' 9i ' Vi) x £ m(Pt2' 92 ■ rft). 

In this equation, q represents the charge of the muon and 
the identities of the two muons are labeled with indices 
1,2. The efficiency of the track selection criteria, e t rk, 
for tracks originating from real muons is determined in 
Ref. using data collected in 2010. Due to the presence 
of additional pp interactions in the same and neighbour- 
ing bunch crossings in 2011 data the systematic uncer- 
tainty on the tracking efficiency has increased, with etrk 
assessed to be 99 ± 1.0% over the whole kinematic range. 

The efficiency to reconstruct a muon, e M , is derived 
using a tag-and-probe method on J/ip — > /i + p~ data. 
The tag muon corresponds to a muon candidate with 
> 4GeV and j^^l < 2.4, and must have fired a singlc- 
muon trigger in the event, as required by the trigger 
matching algorithms. The probe track is only required 
to pass the inner detector track quality, px, and 77 cuts 
and be consistent with having the same vertex as the 
identified tag muon. This technique provides a sample of 
muon candidates unbiased with respect to both trigger 
and offline reconstruction, and with favorable signal to 
background ratio. 

The muon reconstruction efficiency is then derived in 
two-dimensional pr-r] bins (fourteen muon px intervals 
and twenty-six charge-signed muon pseudorapidity [q ■ 77] 
intervals). The ratio of the fitted J/tp —> /i + p~ signal 
yield for those probe tracks identified as muons to the 
fitted signal yield for all probe tracks in these double- 
differential intervals is identified as the single-muon re- 
construction efficiency. 

Due to the toroidal magnetic field, muons with positive 
(negative) charge are bent towards larger (smaller) pseu- 
dorapidity. This introduces a charge dependence in the 
muon reconstruction efficiency. We calculate efficiencies 
as a function of charge-signed pseudorapidity as nega- 
tive muons at positive rapidities arc affected in the same 
manner as positive muons in negative rapidities. The 
charge dependence is particularly noticeable at very large 
| 77 | where the muon can be bent outside of the geometri- 
cal acceptance of the detector. In particular at low px, 
where the muons of particular charge may be bent away 



from (rather than toward for the opposing charge) the 
middle/outer spectrometer planes, they will not be iden- 
tified as muons. 

The efficiency of the di-muon trigger to select events 
that have passed the offline selection criteria, £trig, is also 
calculated from data. This can be factorized into three 
terms: 

Etrig =£Rol(P T l> 9i ' O x £ Roi(p T2 > 92 • ?72 ) 

xc^(AR,\y^\) U 

where £r i is the efficiency of the trigger system to find 
an Rol for a single muon with transverse momentum, 
p T , and charge-signed pseudorapidity, q ■ |t7 m |, and c w is 
a correction for effects related to the di-muon elements 
of the trigger. This correction accounts for the di-muon 
vertex and opposite charge requirements, and for loss of 
efficiency in the di-muon trigger if the two muons are 
close enough together to register only as a single Rol. 

The di-muon correction, c^, itself consists of two com- 
ponents 

c^(AR,y^) = c a (y^) x c AR (AR,y^) (6) 

each evaluated in three separate regions of di-muon ra- 
pidity: barrel (|y MM | < 1.0), barrel-endcap transition 
(1.0 < \y»»\ < 1.2), and endcap (1.2 < \yW\ < 2.25), 
to account for the different behaviors of these corrections 
in these regions. The correction c a is due to the effect 
of vertex and opposite charge requirements on the trig- 
ger and is determined from the efficiency at large di-muon 
angular separations. No difference was observed from de- 
riving the result as a function of A77 and A<f> separately 
rather than versus AR. The asymptotic values are found 
using the ratio of candidate J/ip —> n + n~ decays selected 
by the standard di-muon trigger to those selected by a 
similar di-muon trigger that makes no charge or vertex 
requirements. These values are found to be 99.1 ± 0.4%, 
97.5 ± 0.9%, 95.2 ± 0.4% in the barrel, transition, and 
endcap regions, respectively. 

We extract this dependence on AR in the same three 
regions of di-muon rapidity as used for c Q from a sam- 
ple of offline reconstructed J/-0 — > and Y /i + /x _ 
candidates with p T2 >8 GeV selected using a single-muon 
trigger with a threshold of 18 GeV. The 8 GeV require- 
ment on the lower px muon is made to ensure that the 
efficiency for the trigger system to identify Rols compat- 
ible with muons of px > 4 GeV reached its plateau value. 
The AR dependence of car is then extracted from the 
fraction of fitted di-muon candidates in a AR interval 
within this control sample selected using auxiliary trig- 
gers that additionally pass the di-muon trigger used in 
this analysis. 

The final component of the di-muon trigger efficiency, 
£roI; represents the single-muon trigger efficiency with a 
threshold of p T > 4 GeV. This is measured using well- 
reconstructed J/ip —> fifi candidates in data that pass 
a single-muon trigger (with a threshold of 18 GeV). The 
ratio of the yield of J/ip candidates (determined by fit- 
ting the invariant mass distributions) that pass both this 
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single-muon trigger and the 4 GeV p!^ threshold di-muon 
trigger used in this analysis, to the yield of J/ip can- 
didates that pass the single-muon trigger is identified 
as the single-muon trigger efficiency £r i in the p!^ and 
q- | interval considered. In each case, the reconstructed 
muon(s) are matched to the muon(s) that triggered the 
event for each of the single or di-muon triggers. The num- 
ber of candidates passing the di-muon trigger is then fur- 
ther corrected by c w for di-muon correlation effects [13] . 

The individual and overall efficiencies and weight cor- 
rections, calculated using the methods described above, 
are shown as functions of pt^ and \y^\ in Fig. [3] for the 
T(1S). 



C. Extracting the number of T mesons 

The number of produced Y(nS) mesons used in our 
cross-section determination is found by fitting signal and 
background functions to the m w spectrum of weighted 
Y candidates. The use of per-candidate weights rather 
than average weights allows us to correct for acceptance 
and efficiency without introduction of biases associated 
with the use of average values for these quantities. We 
determine the Y differential cross sections separately for 
each spin- alignment scenario. We perform least squares 
fits to rn^ histograms filled using the weights, w, for 
each candidate in bins of di-muon transverse momentum, 
pt^ 1 , and rapidity, y MM . The form of the \ 2 f° r each (pt^ 1 , 



data 



j=N { 



n data _ n pred 



_2 



(7) 



E 

j=Ni 



where N{ represents the number of T(nS) candidates in 
bin i. 

Predictions for each bin are constructed using four 
sources of di-muon candidates: 



i p ™ d (m^)=N ls F ls (m^;f ls ) 

+ N2sF2s(m PIJ ,; f 2 s) 
+ A r 3S-Fks(m Wi ; / 3 s) 
+ i*b g d(m MM ; /b g d) 



(8) 



where the fx here represent fit function parameters and 
the Fx are normalized probability density functions. In 
order to avoid bias due to the choice of fit model, sev- 
eral parameterizations of signal and background, which 
describe the data well, are used. 

Each of the three Y resonances are parameterized by 
single Gaussian, double Gaussian, or Crystal Ball func- 
tions, chosen to provide a reasonable description of the 
experimental mass resolution and energy loss effects that 
dominate the observed signal line-shapes. The back- 
ground parameterizations vary with di-muon px- At low 
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p^ 1 , and for all rapidity bins, an error function multiplied 
by either a second-order polynomial, or a second-order 
polynomial plus an exponential, is used to model m p)i 
turn-on effects accurately. At mid p 1 ^ a second-order 
polynomial, or second-order polynomial plus an expo- 
nential is adequate to describe the shape of the back- 
ground under the T peaks. At high pt^ a first-order 
polynomial or a first-order polynomial plus an exponen- 
tial is sufficient to describe the background contribution. 
Average fitted values for the numbers of T(nS) mesons 
(N\s, AT 2 g, N$s) taken over all combinations of signal and 
background models are used in the extraction of cross 
sections, while the maximum deviation of any fit result 
from the average is used as an estimate of the fit model 
systematic uncertainty. 

In the fits iVis, A 2 s, A 3 s, the peak mass of the T(1S) 
meson (Mis), parameters describing the shape of the 
Y — > ji + pr mass distribution, and the background pa- 
rameters are allowed to vary freely while M2S, M3S are 
fixed relative to Mis using measured mass splittings [28j j ; 
and the line-shape width parameters <72s, 03s are re- 
lated linearly in to ctis- The distribution of the 
^-probabilities from the fits are uniformly distributed 
indicating that our signal and background parameteri- 
zations provide good descriptions of the data. Figure [4] 
shows the invariant mass distribution of the fitted signal 
region for representative low, mid, and high T p^ inter- 
vals, in both the central and forward rapidity regions. 

In order to account for bin migrations due to finite de- 
tector resolution, corrections in p^f 1 are derived by first 
fitting a smooth analytic function to the background- 
subtracted pt^ 1 spectra of T events in data. This fitted 
distribution is deconvolved with a Gaussian distribution 
using a Gauss-Hermite quadrature integration technique, 
with a pt resolution derived from the fitted invariant 
mass and muon angular resolutions in data. From the 
distributions with and without convolution, migration 
corrections bin- by-bin are derived. These corrections are 
found to be within 0.5% of unity for central rapidities 
and within 1% of unity for forward rapidities. The num- 
ber of efficiency-corrected T —> /i + fi~ decays extracted 
from the fits in each (pjf 1 , bin are corrected for the 
difference between true and reconstructed values of the 
di-muon and rapidity. Using a similar technique, bin- 
migration corrections as functions of are found to be 
negligible. 



V. SYSTEMATIC UNCERTAINTIES 

We consider the following sources of systematic uncer- 
tainty on the Y(nS) differential cross sections: luminos- 
ity determination, reconstruction and trigger efficiencies, 
migration between pt^f and bins due to resolution, 
acceptance corrections, and the background and signal 
fit models used. The range of these uncertainties for 
the three Y states is summarized in Table U and their 
breakdown in each source is given for the corrected cross- 



section analysis in Figs. EH7J 

The relative luminosity uncertainty of 3.9% is de- 
scribed in more detail in Ref. [29j . Other sources of 
systematic uncertainty are discussed below. As the sta- 
tistical components of the uncertainties associated with 
the determination of e rcco and £t r i g are dominant, the 
uncertainties on the cross sections are derived from the 
statistics of the control samples used to extract them us- 
ing a series of pseudo-experiments, randomly varying the 
weights used for each Y candidate in data, according to 
the uncertainties on the efficiency maps. Systematic un- 
certainties associated with the fit model used to extract 
the number of Y(nS) decays from our di-muon data sam- 
ple, are quantified by taking the largest deviation of the 
fitted values of N n s found in the six possible model com- 
binations from the average value. Systematic effects, due 
to differences in the underlying kinematic distributions of 
the control sample and the data distributions, are found 
to be negligible due to the fine differential binning of the 
efficiency and acceptance corrections. 

Uncertainties associated with the acceptance correc- 
tion include statistical uncertainties on the determina- 
tion of the correction in fine pt and rapidity bins. This 
constitutes approximately a 0.5% uncertainty across the 
measured spectrum. A shift in the interaction point along 
the beamline axis, z, influences the acceptance particu- 
larly at large rapidity. To estimate the impact of this ef- 
fect, the acceptances were recalculated with shifts in z by 
±62 mm, corresponding to one Gaussian standard devia- 
tion of the vertex z position distribution in the analyzed 
data, and the variation in the measured cross sections 
was calculated as a function of Y p^ and rapidity. These 
variations result in uncertainties on the acceptance cor- 
rections of 0.4-0.7% as a function of pt, and 0.6-1.5% 
as a function of rapidity, growing toward larger rapidi- 
ties and lower transverse momenta. Our estimate of the 
correction to the number of Y — > decays fitted in 

each (p!p\ bin for differences between reconstructed 
and true values of p 1 ^ due to bin migration depends on 
a knowledge of the di-muon px resolution. Allowing this 
resolution to vary within its uncertainty results in a neg- 
ligible change in the Y(nS) cross sections. 

As mentioned earlier, Y spin-alignment effects have an 
impact on the determination of cross sections extrapo- 
lated to full phase space. As in the ATLAS measurement 
of J/tp production [7J , we quote an uncertainty due to Y 
spin-alignment by comparing results obtained using the 
unpolarized assumption (Ag = A^ = Ag0 = 0) to those 
derived under the other four extreme parameter combi- 
nations. The size of the possible variation in each bin 
is taken to be the largest positive and negative devia- 
tion from the unpolarized baseline result. We provide 
the cross-section results for each of the spin-alignment 
scenarios along with tabulated values [26| of the weight 
corrections so that the corrected data can be determined 
in terms of any spin-alignment scenario. 
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TABLE I: Summary of statistical and systematic uncertainties on the cross-section measurements. The systematic uncertainty 
due to the acceptance determination applies only to the corrected cross-section measurements and does not include possible 
variation of the result due to spin-alignment. Values quoted refer to the range of uncertainties over in each of the 
regions considered. 



Relative Uncertainty [T(1S) / T(2S) / T(3S)] [%] 

Source < 1.2 1.2 < \y^\ < 2.25 

Statistical Uncertainty 1.5-10 / 2-10 / 3-10 2-15 / 2.5-12 / 4-18 

Fit Model 0.3-2 / 0.3-10 / 0.3-15 0.4-10 / 0.3-5.0 / 0.2-12 

Luminosity 3.9 3.9 

Trigger Efficiency 2.0-3.0 / 1.0-2.0 / 1.0-2.0 2.0-3.0 / 2.0-4.0 / 2.0-7.0 

Muon Reconstruction Efficiency 0.7-2.0 / 0.4-1.2 / 0.4-1.2 0.5-5.0 / 0.4-3.0 / 0.4-3.0 

Acceptance corrections 0.7-1.5 0.7-1.5 

Track Reconstruction Efficiency 1.0 1.0 
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central rapidity region. A common luminosity uncertainty of ±3.9% is not included. 
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VI. RESULTS AND DISCUSSION 



Fiducial cross sections 



We measure the differential cross section multiplied by 
the di-muon branching fractions of T(1S), T(2S), and 
T(3S) mesons as a function of T transverse momen- 
tum and rapidity, both in a fiducial region defined by 
pt^ > 4 GeV, \rj^\ < 2.3 (free from spin-alignment uncer- 
tainties), and corrected back to the full muon decay phase 
space with \y r \ < 2.25. We additionally present mea- 
surements of the production cross sections of the T(2S) 
and T(3S) relative to the T(1S) as a function of T(1S) 
transverse momentum and rapidity. Tabulated results of 
all measurements presented in this paper are available in 

HEPDATA [26| . 



Differential cross sections multiplied by the T fi + fj," 
branching fractions, d 2 a / 'dp^dy x Br(T — > p + n~), are 
calculated within the fiducial acceptance of our analy- 
sis (p£ > 4 GeV, \if\ < 2.3) using the results of fits 
to di-muon candidates in data corrected, candidate by 
candidate, for efficiencies using the fiducial weights, tugd, 
defined in Eq. ((2J. These differential fiducial cross sec- 
tions along with total uncertainties are shown in Figs. [5] 
and |9] and span four orders of magnitude across the pt 
range studied. 

Integrating over all p\ and both rapidity bins we find 
cross sections within the fiducial acceptance of the detec- 
tor as shown in Table |TTJ The results for the T(1S) are 
consistent with our previous measurement [8[ but extend 
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FIG. 7: Sources of statistical and systematic uncertainty on the corrected T(nS) production cross-section measurements as a 
function of T rapidity. A common luminosity uncertainty of ±3.9% is not included. 



to a significantly higher p^ and with increased precision. 



TABLE II: Integrated fiducial cross-section measurements for 
T(nS). Uncertainties quoted represent statistical, systematic, 
and luminosity terms respectively. 

Integrated fiducial cross sections 

p£ >4GeV, |^| < 2.3 

State crfld(pp — > T) x Br(T — > n + n~) 

Range: pi < 70 GeV, \y r \ < 2.25 

T(1S) 1.890 ± 0.007 ± 0.095 ± 0.074 nb 
T(2S) 0.601 ± 0.003 ± 0.040 ± 0.023 nb 
T(3S) 0.304 ± 0.003 ± 0.021 ± 0.012 nb 



Presenting the results in a restricted kinematic phase 
space removes any uncertainty due to the spin-alignment 
of quarkonia. This allows unambiguous comparison 
of differential spectra with any theoretical approaches 
which can provide predictions with kinematic restrictions 
applied to the decay products of the T. 

NLO color-singlet [14] predictions [3(| have previously 



been compared to differential T(1S) fiducial production 
cross sections [1[ and were found to underestimate the 
measured production rates by approximately an order of 
magnitude and to not reflect the pt dependence of the 
data. This is not surprising as it has been known [l6| for 
some time that higher-order corrections are both large 
and necessary in order to adequately describe quarko- 
nium production at high with color-singlet terms. 
At this time, no color-singlet calculations beyond NLO 



12 



> 

CD 

O 

Q. 



T 

m 

x 
>, 

T3 

h 
Q. 
"O 

e 

CM 

T3 



10* 
10 

1 

.10" 1 
10" 2 



' 1 I 1 11 1 I 1 11 1 I 11 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 

Fiducial cross sections (p|* > 4 GeV, < 2.3) - 

|y r |<1.2 

gjjjYCIS) 
M T(2S) 

U V T(3S) 



■A- • 



A A>. 



y<s = 7TeV j Left =1.8 fb" 1 



' 10 " 20 " ' 30 ' 40 50 60 70 



Tp T [GeV] 



ID 
O 
O 
23 
Q. 

+ 

T 

m 

x 

~B 

T3 



10 c 



10' 



10 



_ , , , , I , , , , I , , 


I I - 


_ Fiducial cross sections (p 


> 4 GeV, |rf| < 2.3)_ 




" * 'Mm 




* ■ 
i 


ATI AG 

V r s = 7TeV f Ldt = 1.8 fb 1 


: p{: < 70 GeV 




I »T(1S) 




I T(2S) 




T(3S) 




: , , , , 1 , , , , 1 , , 


,i,ii,i,,= 



0.5 



1.5 



T|y| 



> 

CD 
D. 



*2 10 

T 

£ 1 

CD 

X 



Fiducial cross sections (p T > 4 GeV, |rfj < 2.3) - 

K1.2<|y r |<2.25 
■ rtis, 

|T(2S) 

■ T(3S) 



■D H 10 
Q. 

b io _ - 



I , , i i | 



* * • 



V 



l=s = 7TeV J L dt = 1 .8 fb 



10 20 30 40 50 60 70 

Tp T [GeV] 

FIG. 8: Differential cross sections multiplied by the di-muon 
branching fraction, d 2 a /dpTdy x Br(T —¥ for cen- 

tral (top) and forward (bottom) rapidities, for T(1S), T(2S), 
T(3S) production within the fiducial acceptance. Points with 
error bars indicate results of the measurements with statisti- 
cal uncertainties while shaded areas correspond to total un- 
certainties on the measurement, including systematic effects. 



are available for quarkonium production measurements 
quoted in a restricted muon acceptance, nor are Color 
Octet [13] or Color Evaporation [l9j approaches currently 
able to account for the kinematics of leptons from the 
quarkonium decay. The fiducial measurements presented 
here are precise and free from any assumptions on the 
angular dependencies of the di-muon system, accurately 
reflecting the production dynamics in proton-proton col- 
lisons at ^/s = 7 TeV. 



FIG. 9: Differential cross sections multiplied by the di-muon 
branching fraction, da/dy x Br(T — > /i + /i~)(pT-integrated) 
da/dy x Br(T -> n + /i,-)far T(1S), T(2S), T(3S) production 
within the fiducial acceptance. Points with error bars indi- 
cate results of the measurements with statistical uncertain- 
ties while shaded areas correspond to total uncertainties on 
the measurement, including systematic effects. 



B. Corrected cross sections 

We also calculate differential cross sections multiplied 
by the T — > branching fractions extrapolated to 

the full muon phase space within \y r \ < 2.25. For these 
corrected cross sections, the results of fits to di-muon 
candidates in data are corrected, candidate by candi- 
date, for efficiencies using the total weights, w to t, de- 
fined in Eq. © . Results are shown for the isotropic spin- 
alignment scenario in Figs. 1101 and [Til as a function of T 
Pt and y for all three states and, integrated over all pi^ 
and both rapidity bins, in Table [TTT1 

Our results are consistent with measurements by 
CMS @ and, in the small region of rapidity overlap from 
2.0 < \y\ < 2.25, with LHCbflU- These results allow 
us to test phenomenological models of Y production not 
just at low pt (for complementarity with results from the 
Tevatron experiments) but in a newly-probed region of 
significantly boosted T, where higher-order contributions 
become particularly important. 

Figures [T2"rfT7l show the differential cross sections as 
a function of pt and rapidity for each of the T states 
in comparison with theoretical predictions. The effects 
of varying spin-alignment assumptions from the nominal 
assumption of isotropic muon angular distributions inde- 
pendent of 6* and (f>* are indicated by a shaded band. 

Clearly spin- alignment has a large effect on the T pro- 
duction cross sections, especially at lowp^, and in partic- 
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FIG. 10: Differential cross sections multiplied by the di-muon 
branching fraction, for T(1S), T(2S) and T(3S) production 
extrapolated to the full phase space for the (top) \y T \ < 1.2, 
(bottom) 1.2 < \y r \ < 2.25 rapidity intervals. Points with 
error bars indicate the results of the measurements with total 
statistical and systematic errors. Results are shown assuming 
an isotropic spin-alignment scenario. 



ular if there is a non-trivial azimuthal component to the 
spin- alignment. New results [3, [l2| of T spin-alignment 
from CDF and CMS suggest that the spin-alignment is 
consistent with unpolarized production. Our central as- 
sumption of isotropic T decays is consistent with these 
results. Nevertheless, as these spin-alignment measure- 
ments are made at different center-of-mass energies or 
in a restricted phase space in both p^. and rapidity with 
respect to measurements presented here, we provide the 



FIG. 11: Differential cross sections multiplied by the di-muon 
branching fraction, for T(1S), T(2S) and T(3S) production 
extrapolated to the full phase space, p-r-integrated as a func- 
tion of absolute T rapidity. Points with error bars indicate 
the results of the measurements with total statistical and sys- 
tematic errors. Results are shown assuming an isotropic spin- 
alignment scenario. 



TABLE III: Corrected cross-section measurements in the 
isotropic spin-alignment scenario. Uncertainties quoted rep- 
resent statistical, systematic, and luminosity terms, respec- 
tively. 

Integrated corrected cross sections 
State a(pp -> T) x Br(T -> a* + /0 

Range: p£ < 70 GeV, \y r \ < 2.25 

T(1S) 8.01 ±0.02 ±0.36 ±0.31 nb 

T(2S) 2. 05 ±0.01 ±0.12 ±0.08 nb 

T(3S) 0.92 ± 0.01 ± 0.07 ± 0.04 nb 



results under a variety of polarization scenarios so that 
the impact of spin- alignment on the corrected cross sec- 
tions can be quantified across the full range of study. 

The contributions of the five polarization scenarios can 
be seen in the lower panes of each plot where the ratio of 
the differential cross section under these spin-alignment 
assumptions to the unpolarized scenario is shown. Across 
the whole pr range studied the envelope is bounded from 
above by the T±± (Xg = + l,X$ = ±1, Xg^ = 0) scenario 
with a maximal </>* variation. From below, the cross- 
section envelope is bounded by fully longitudinal spin- 
alignment at very low pt, with the T-\ — (Xg = ±1, A^ = 
— l,Ae0 = 0) scenario resulting in the largest downward 
variation at pt > 4 GeV. In this measurement we extend 
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FIG. 12: Differential cross sections multiplied by the di-muon 
branching fraction, d 2 a / dp^dy x Br(T — s> /i + /i~), for T(1S) 
production extrapolated to the full phase space for (top) cen- 
tral and (bottom) forward rapidities. Points with error bars 
indicate results of the measurements with total (statistical 
and systematic) uncertainties. The maximal envelope of vari- 
ation of the result due to spin-alignment uncertainty is indi- 
cated by the solid band. Also shown are predictions of di- 
rect production with the NNLO* Color Singlet Mechanism 
(CSM) and inclusive predictions from the Color Evaporation 
Model (CEM). Direct CSM predictions are corrected to inclu- 
sive production predictions by a multiplicative factor of 1.5 
to account for the effects of feed-down. These theory predic- 
tions are shown as a ratio to the data in the lower panes for 
CEM (middle) and CSM (bottom), along with detail of the 
variations of the cross-section measurement under the four 
anisotropic spin-alignment scenarios as a ratio to the nominal 
data. 



FIG. 13: As for Fig.[l2l but for T(2S) production. No correc- 
tion is made to the direct CSM predictions for the unknown 
feed-down contribution in this case. 



the pt range above 45 GeV where the maximal possible 
impact due to the unknown spin-alignment of T is below 
±10%. This is significantly smaller than the theoreti- 
cal uncertainties and is of similar magnitude to current 
experimental uncertainties. As such, this region will of- 
fer a precision environment to compare future theoretical 
studies of T production to data. 

In Figs. [12]-[T4] a comparison is also made to two the- 
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FIG. 14: As for Fig.[l2l but for T(3S) production. No correc- 
tion is made to the direct CSM predictions for the unknown 
feed-down contribution in this case. 



oretical predictions of Y production. The first [16| , is 
a QCD-based calculation using the color-singlet mech- 
anism referred to as NNLO* CSM and presumes 
that T meson production occurs via a color-singlet state 
and includes full corrections up to next-to-leading order 
(NLO), as well as some of the most important next- 
to-next-to-leading-order (NNLO) terms. This inclusion 
significantly modifies the prediction. The partial na- 



ture of the higher-order calculation limits the applicabil- 
ity of the calculation to values above a particular T pt 
threshold and increases the sensitivity of the prediction 
to the choice of renormalization and factorization scales. 
The second prediction, known as the Color Evaporation 
Model [_01 Hlj], labeled as CEM, is a phenomenological 
model for inclusive Y production based on quark-hadron 
duality. This model assumes that any heavy QQ pair pro- 
duced in the initial hard scattering evolves to a quarko- 
nium state if its mass is below the threshold of open 
heavy-flavor meson pairs. Predictions of the CEM model 
involve a single constant that must be determined from 
corrected cross-section measurements for each quarko- 
nium state and use a 6-quark mass of 4.75 GeV. As in 
the case of the CSM, divergences in the predictions re- 
strict the applicability of the model at low Y p^. Uncer- 
tainties from factorization and renormalization scales are 
estimated from varying the scales independently up and 
down by a factor of two and additional uncertainties are 
estimated from varying the 6-quark mass [HJ ■ 

Color-singlet calculations at NLO/NNLO* predict a 
largely longitudinal polarization of direct Y particularly 
at high pr (although the effect of feed-down introduces 
large uncertainties to this prediction). The Color Evap- 
oration Model offers no explicit prediction of the spin- 
alignment evolution of Y but the nature of the model 
suggests no strong polarization should be observed as no 
single production mechanism dominates. 

As can be seen in Figs. [T^] HU the two models provide 
quite different descriptions of Y production. Predictions 
from CSM are for direct Y production only, and so do not 
account for Y production that arises from feed-down from 
the production of higher Y states or from radiative decays 
of the Xf,j(nP). From a previous measurement (33j the 
contribution of feed-down to T(1S) production is known 
to be approximately 50%. We correct the CSM T(1S) 
prediction by scaling the direct prediction by a factor of 
1.5 (no such correction is needed for the CEM as this is 
already an inclusive calculation). 

Feed-down contributions are expected to similarly con- 
tribute to T(2S) production, but no measurement or 
reliable prediction for the relative contribution exists, 
so we do not apply a correction to the direct T(2S) 
CSM predictions. For T(3S) there is no feed-down from 
higher Y states. Recently, the ATLAS experiment discov- 
ered the existence of a state or states interpreted as 
the Xbj(3P) below the BB threshold that are expected to 
have a significant branching fraction for radiative decays 
into T(3S) + 7 and thus induce a feed-down contribution 
to T(3S) (and other Y states). As the relative production 
and decay rates of these states are also as yet unknown 
no correction is applied to the CSM predictions for T(3S) 
production either. 

For each of the three Y(nS) states, the NNLO* CSM 
predictions (considering also the additional normaliza- 
tion uncertainty due to the poorly-known contributions 
from feed-down) fit our data well in the moderate p\ 
region, but exhibit a steeper px dependence than seen 
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in data. The predictions therefore significantly under- 
estimate the cross section at high pr- Theoretical de- 
velopments in the prediction of feed-down contributions 
may improve this description. CEM predictions appear 
to show a better match with data at high These 
predictions underestimate the rate (favouring a smaller 
choice of renormalization/factorization scale or larger b- 
quark mass) and have problems in modeling the shape of 
the spectrum, particularly at lower px- Shape discrep- 
ancies cannot be accommodated within the uncertainties 
quoted as changes in the scale choice introduce correlated 
changes in the prediction as a function of pr ■ 
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FIG. 15: Differential cross sections multiplied by the di-muon 
branching fraction, da/dy x Br(T — > for T(1S) pro- 

duction extrapolated to the full phase space. Points with error 
bars indicate results of the measurements with total (statis- 
tical and systematic) uncertainties. The maximal envelope of 
variation of the result due to spin-alignment uncertainty is in- 
dicated by the solid band. The variations of the cross-section 
measurement under the four anisotropic spin-alignment sce- 
narios are shown in the lower pane as ratios to the unpolarized 
scenario. 
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FIG. 16: As for Fig. [I5j but for T(2S) production. 
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Figures [15] to [17] show the variation of the production 
cross section as a function of absolute T rapidity inte- 
grated across all px- The dependence on rapidity is rel- 
atively flat in the interval of rapidities in the ATLAS 
acceptance. 

Variations between spin-alignment scenarios lead 
largely to a change in the normalization of the distri- 
butions with little variation in shape except at high ra- 
pidity. There, the T-\ — I- scenario leads to an increase in 
the cross section with increasing rapidity, while the fully 
longitudinal scenario leads to a drop in cross section at 
high (|y| > 1.7) rapidity. 



FIG. 17: As for Fig. [HI but for T(3S) production. 



C. Cross section ratios 

From our differential cross-section measurements, we 
explore the pr and rapidity dependence of T(2S) and 
T(3S) production relative to the T(1S) by deriving the 
ratios 

d 2 a(pp -> T(nS))/d PT dy ^ Br(T(nS) -» {SjjT) 
d 2 <r(pp -S- Y{lS))/d PT dy * Br(T(lS) ->■ ( ' 
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for \y | < 2.25 with n=2,3. Such observables are sen- 
sitive to the magnitude and kinematic dependencies of 
feed-down contributions between the three T states. Re- 
sults of the differential cross-section ratio measurements 
are presented in Figs. [T51and[T51 
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FIG. 18: Ratios of differential T(2S)/T (IS) and T(3S)/T(1S) 
cross sections multiplied by the di-muon branching fractions 
as a function of di-muon rapidity. Points with error bars in- 
dicate results of the measurements with statistical uncertain- 
ties while shaded areas correspond to total uncertainties on 
each point, including systematic effects, but excluding spin- 
alignment effects. 



These measurements are made under the assumption 
of unpolarized T(nS) mesons, and take into account sta- 
tistical correlations between the fitted numbers of T(1S), 
T(2S) and T(3S) mesons. Systematic uncertainties are 
estimated by varying acceptance, efficiency, and fit model 
assumptions coherently in the numerator and denomina- 
tor when calculating the ratios, thereby partially cancel- 
ing uncertainties in the ratio. Luminosity uncertainties 
cancel entirely. 

The measured T(nS)/T(lS) ratios are relatively con- 
stant in the < pt < 5 GeV interval at ~ 20% and 
~ 7% respectively for the T(2S) and T(3S). At higher 
Pt a significant and steady rise in the relative produc- 
tion rates of higher T states is apparent, in agreement 
with measurements by CMS@. However, at the larger 
Pt values (above p^ of 30-40 GeV) accessible for the first 
time with these measurements, evidence of a saturation 
in this rise is apparent, suggesting that we are probing a 
regime where direct production dominates over contribu- 
tions from the decays of excited states. In contrast, the 
rapidity dependence of these production ratios is quite 
flat across the full \y\ < 2.25 rapidity interval. 

Higher-order color-singlet calculations are not cur- 
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pidity regions. Points with error bars indicate results of the 
measurements with statistical uncertainties while shaded ar- 
eas correspond to total uncertainties on each point, including 
systematic effects, but excluding spin-alignment effects. 



rently able to predict the evolution of these production 
ratios in p^, due to the significant feed-down contribu- 
tions. At leading order in the quark velocity in the per- 
turbative expansion, the production ratio of the direct 
contributions is proportional to the ratio of the squares 
of the magnitudes of the wavefunction at the origin (or 
the partial decay widths), multiplied by the branching 
fractions to di-muons for each of the states in question. 
The predicted ratios then are 36% for T(2S)/T(1S) and 
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29% for T(3S)/T(1S) for direct production. At high p T , 
where the measured production ratio plateaus, the val- 
ues are somewhat higher than these predictions suggest- 
ing that there is some enhancement (with respect to the 
simple picture above) of the T(2S) and T(3S) production 
relative to T(1S) production. 

VII. CONCLUSIONS 

We have measured differential production cross sec- 
tions and relative production rates for T(1S), T(2S), 
and T(3S) mesons in pp collisions at y/s = 7TeV at 
the LHC up to p^. < 70 GeV in the rapidity inter- 
val \y r \ < 2.25. The possible impact of the T spin- 
alignment on these measured spectra is also quantified. 
These measurements are compatible with measurements 
by the CMS and LHCb collaborations. We have mea- 
sured integrated corrected cross sections multiplied by 
the T — >• /i + p~ branching fractions, within the rapid- 
ity region \y r \ < 2.25 of 8.01 ± 0.02 ± 0.36 ± 0.31 nb, 
2.05±0.01±0.12±0.08nb, and 0.92±0.01±0.07±0.04 nb 
for the T(1S), T(2S), and T(3S), respectively. Uncertain- 
ties correspond to statistical, systematic, and luminosity 
measurement effects. 

These cross-section estimates assume unpolarized pro- 
duction corresponding to Xg = X^ = Xg^ = in 
Eq. ([3]). If the production polarization is fully trans- 
verse or longitudinal with no dependence on the az- 
imuthal decay angle 0*, corresponding to Xg = ±1, = 
Ae,/, = 0, the integrated cross sections may vary by up 
to (+19,-23)%, (+18,-21)% and (+17, -19)% respec- 
tively for the T(1S), T(2S) and T(3S). However, accep- 
tances corresponding to certain non-zero values of X^ can 
become extremely low, especially at low transverse mo- 
menta of the muon pair. This would lead to a significant 
increase of the corrected cross sections in these configu- 
rations, up to +217%, +155% and +126% for the three 
T states. 

We compare our differential cross-section results to 
predictions from two theoretical approaches describing 
T production. Our measurements find both the NNLO* 
CSM and the CEM predictions have some problems in 
describing the normalization and shape of the differen- 
tial spectra. In particular, NNLO* dramatically under- 
estimates the rate at high transverse momenta, where 
the data tend to agree better with the CEM. The in- 
clusion of P-wave feed-down contributions in the theo- 
retical calculation may help to improve the description. 
Large scale uncertainties in these predictions allow pos- 
sible contributions from color-octet terms to contribute 



to the production rate in addition to singlet diagrams. 
The differential production ratios indicate that the in- 
crease in the production of higher T states as a function 
of relative to the T(1S) observed previously begins to 
saturate at 30-40 GeV. Above ~ 40 GeV, the envelope of 
possible variations in the differential cross sections due 
to spin-alignment is reduced to below ±10%. This, along 
with the expected reduction in feed-down contributions, 
results in a relatively well-controlled region in which to 
study quarkonium production without the dominant ex- 
perimental and theoretical effects that complicate such 
studies at lower p^. 
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